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ABSTRACT 



Many scientific and technical endeavors require the reconstruct ion 
of a three-d irnensional solid from a collection of two-d irnensional 
contours. One method for this reconst ruct ion involves a procedure 
whereby individual pairs of contours are mapped together to form 
triangular surface patches. In this paper, we present an algorithm 
which not only handles mapping situations of simple, closed contours 
but also mappings of multiple contours per plane and partial contour 
mappings. Also included is a discussion of algorithm limitations and 
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I. INTRODUCTION 



Nan y scientific and technical endeavors require the reconstruct ion 
of a three-dimensional solid from a collection of two-dimensional 
planar contours. These contours are obtained by sorne sensor method 
that samples the original three-dimensional solid along a finite number 
of parallel planes. The data extracted from that set of parallel 

planes are contours that lie along the solid’s exterior and interior 
surfaces. The contours on the parallel planes appear as line segments. 
The line segments are either closed loops, open segments, or single 
points. The goal of surface construction is the formation of surface 
patches between contours on adjacent planes such that an approxirnat ion 
of the original three-dimensional solid is formed. 

Surface construct ion by the tr iangulat ion of two-dimensional 
contours is the procedure by which a pair of parallel, planar contours 
are "mapped together" and then "triangulated" into surface patches that 
form a surface display. The mapping operation of the surface 
construction algorithm identifies which contours on consecutive, 
parallel planes should be mapped together, and exactly which portions 
of those contours should be connected. The tr iangulat ion operation 
forms the connections between contours on adjacent planes by building 
triangular tiles between those mapped contours. Each triangular tile 
is built from an individual line segment from one contour and a single 
point from the end of a line segment on the other mapped contour. This 
tiling operation is performed for all line segments in the connect 
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region of each mapped contour. The connect region is that section 



of coordinates designated as mappable for a pair of contours on 
consecutive planes. 

Not at ional ly, this problem has been specified as follows: 

"An unknown three dimensional solid is intersected by a finite 
number of specified parallel planes. 

The only information about the solid consists of the 
intersect ions of its surface with the planes. Each of these 
intersect ions is assumed to be a simple closed curve. These curves 
are not completely specified; instead, a finite sequence of points 
encountered during a positive (counterclockwise) traversal of each of 
the original curves is given. The curve segment between two 
consecutive points is approximated by a linear segment, called a 
contour segment. 

We reduce the problem of constructing such an approximating 
surface to one of constructing a sequence of partial approximat ions, 
each of them connecting two contours lying on consecutive planes. 
[Figure 1.1] 

Let one contour be defined by the sequence of rn distinct contour 
points P0, PI, ..., P(m-l), and let the other contour be defined by 
the sequence of n distinct contour points Q0, Ql, ..., Q(n-l). We 
note that P0 follows P(rn-l) and that Q0 follows Q(n-l), and so 
indicies of P are modulo rn and indicies of Q are modulo n. We wish 
to create a surface between the contours P and Q. The surface is 
constructed of triangular tiles between these two contours. The 
verticies of these tiles are contour points, with the verticies of 
each tile taken two from one sequence and one from the other. Thus, 
each tile is defined by a set of three distinct elements either of 
the form -CPi,Pk, Qj> or -CG!i,Qk, Pj>. [Figure 1.23 

Each tile’s boundary will consist of a single contour segment and 
two spans, each connecting an end of the contour segment with a 
common point on the other contour. " [Ref. 13 

This not at ional specification of the problem is consistent in all 
papers accessible in the literature on surface construction [Ref. 13 
[Fef . £3 [Ref. 33 [Ref. 43. 
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Fig 1.2 - Mapped connections into triangulated 
surf ace patches . 
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The initial emphasis of this paper is a review of the previous 
algorithms for surface construct ion. Included in this review is a 
discussion of each algorithm’s capabilities and limitations. After 
this review, we present a new algorithm for surface construction that 
is more comprehensi ve than any that has previously appeared in the 
literature. Following that discussion, we examine the limitations of 
our new algorithm. 
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II. LITERATURE REVIEW 



In order to understand the nature of our new algorithm for surface 
construction, it is important that we understand the previous 
algorithms for surface construct ion. Four such algorithms have 
provided the background necessary for the development of our algorithm. 
These are presented in chronological order. 

A. FUCHS ALGORITHM 

The first algorithm we examine for the reconstruct ion of a 
three-dimensional object from its planar contours was presented by 
Fuchs in 1377 CRef, 13. The problem statement from that article 

(reproduced in our introduct ion) has been used in all subsequent papers 
which build upon the Fuchs algorithm. The major contribution of that 
article, in addition to the concise statement of the problem, is the 
present at i on of an algorithm capable of connecting simple, closed 
contours (Figure 2.1). 

The problem with the Fuchs algorithm stems from its inability to 
handle multiple contours on adjacent planes (Figure £.2). 

Add 1 1 ional ly, no mechanism is provided to handle oartial contour 

mappings or open (non-closed) contours. With respect to the case of 

multiple contours on adjacent planes, no mechanism is provided to 

identify which of the contours should be mapped together. The general 

case for surface construct ion is to have multiple contours on each 
plane. The problem with partial contour mappings is that the Fuchs 
algorithm can only construct a complete tr iangulat ion between adjacent 
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Fig. 2.2 - Example of multiple contours per plane. 
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contours. 



This limitation disallows partial tr iangulat ions of 



contours. Such partial mappings often are indicated for cases of 
dissimilarly sized contours. Finally, the problem of open contours can 
be attributed to algorithm generality. A mechanism that solves the 
partial contour mapping problem can also solve this problem. 

B. CHRISTIANSEN ALGORITHM 

In the Christiansen paper, an algorithm is presented which is 
similar to the Fuchs algorithm. The major dissimilarity is the 
inclusion of a mechanism to facilitate human interaction for the 
resolution of highly ambiguous contour mappings. Human interaction is 
used to determine the relative connection points in the contour mapping 
process for highly convoluted contours. 

Similarly to the Fuchs algorithm, this algorithm can handle 
mappings of simple, closed contours. It also has capabilities for 
mapping together simple branches. An example of such branching, seen 
in Figure 2. 3, is a pair of contours on one plane being mapped to a 
single contour on an adjacent plane. This capability allows the 
algorithm to handle simple cases of multiple contours on adjacent 
planes. The method by which this problem is solved is as follows: 

1. Introduce a new node midway between the closest nodes on the 
branches. The Z coordinate of this node is the average of the Z 
coordinates of the two contour levels (planes) involved. 

2. Renumber the nodes of the branches and the new nodes such that 
they can be considered as being one loop. [Figure 2. 4] 

3. Triangulate as usual. [Ref. 2: pp. 189-190] 
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Fig. 2.4 - Triangulation scheme for branching. 
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The Christ i arisen algorithm is not capable of handling open contours, 
nor is it capable of handling complex cases of multiple contours on 
adjacent planes, except by way of expensive human interaction. ft final 
note of interest with respect to this algorithm, is the use of a 
heuristic for selection of the nodal connections. In cases where 
contours on adjacent planes are mutually centered and are reasonably 
similar in size and shape, selection for nodal connection is based on 
"shortest diagonal" rather than minimum triangular area [Ref. 2: p. 
1883. During this operation, one of two nodes is selected to create 
the next triangular surface patch. The nodes under considerat ion are 
the two "next" nodes of each contour. 3y determining the length of 
each of the possible diagonals for the surface patch, the connection 
node is selected based on minimum length. 

C, SnftNTZ ftLSORITHM 

The algorithm presented in the Shantz article CRef. 33 extends the 
algorithms of Fuchs and Christiansen to handle contour defined objects 
which are highly branched and have holes. Multiple contours on 
adjacent planes are handled by 

", . . first concatenat ing the contours on each plane into a 

single large contour using minimum distance links, then performing 
the mapping between the resulting composite contours. " [Ref. 3: p. 
2423 

Shantz uses the simple, closed contour mechanism of Fuchs to form tne 
connections between the composite contours. Once the connections have 
been formed, the extraneous ones (due to concatenat ion) are removed. 
Some difficult multiple contour cases for this algorithm require human 
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interaction to 



solve ambiguities. 



Similar to the Christiansen 



algorithm, Shantz states that this should be avoided since human 
interaction is "extremely labor intensive. " He cites a case which 
required 50 to Q0 hours of contour splitting, using an interactive 
cursor, to produce a surface display for the highly convoluted cortex 
and basal ganglia contours (extracted from the Livingston brain 
dat abase) . 

This algorithm, similar to the Christiansen algorithm, is limited 
in its ability to handle cases of open contours and partial contour 
mappings. PI so, cases of multiple contours on adjacent planes can be 
handled only when a composite contour can be formed, or when 
ambiguities are resolved via human interaction. 



D. GPNPPPTHY ALGORITHM 

The most recent algorithm for surface construction from planar 
contours was presented in a paper by Ganapathy CRef. 4]. That 
algorithm is essentially an improvement on the Fuchs and Christiansen 
algorithms for simple, closed contours, without the capabilities 
described by Shantz. Like Fuchs, Ganapathy assumes a complete mapping 
of contours, which is not always possible. The improvement over the 
Fuchs and Christiansen algorithms is attributed to the use of a more 
computationally expedient heuristic for triangulat ions. 

The problem with the Ganapathy algorithm is that it presents a 
general solution for handling only the simple case of mapping single, 
closed contours on adjacent planes. The issues of multiple contour 
mappings and partial contour mappings are ignored. Additionally, no 



18 



mechanism for user interaction is provided for resolving mapping 
ambiguities, further limiting the algorithm to simple cases. 

E. SUMMARY 

None of the above papers provides a complete solution to the 
problem of surface construction via the tr iangulat ion of contours. 
What is required is an algorithm with capabilities for multiple 
contours per plane and partial contour mappings. Additionally, the 
algorithm should support simple cases of branching and provide a 
mechanism for human interaction for the resolution of highly ambiguous 
mappings. 

The surface construction algorithm we present handles not only the 
simple contour mapping problem, but also provides a more comprehensive 
procedure for solving the multiple contours per plane and partial 
mapping problems. The only capability lacking from our algorithm is 
that for handling branching as per the Christiansen paper. A detailed 
discussion of our algorithm follows, with a presentation of algorithm 
heuristics and limitations. 
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III. THE ALGORITHM 



In the last section, we presented a discussion of previous 
algorithms for surface construction via triangulation. Here, we 
present an in-depth discussion of our algorithm by first discussing 
known input /out put data structures. Following this presentat ion, an 
overview of the major parts of the algorithm precedes a detailed 
discussion of the parts. 



A. INPUT /OUTPUT SPECIFICATIONS 

The problem of surface construction of an object from a set of 
planar contours, as seen in Figure 3.1, can be reduced to one of 
constructing the surface tr iangulat ions between two adjacent planes. 
The specification of the problem can be best viewed by detailing the 
known input data structures: 



* t ot a 1 ( i ) 

* start ( j, i ) 

* length ( j , i ) 

* type ( j , i) 

* interior ( j , i ) 



number of contours on plane i. 
start of contour j on plane i. 

number of coordinates in contour j on plane i. 

type of contour j on plane i. 

(CLQSED_LQOP, QPEN_SEGMENT, or SINGLE_POINT) 

value of contour j’s interior with respect to 
the contour line. 

(HIGH, LOW, or INDETERMINATE) 



+ coords (XYZ, pointer, i ) : input coordinates for all contours on plane 

i. To isolate contour j on plane i: 

for (pointer = start (j, i) + k - 1), 
where k = 1, length(j,i). 



From the above data, we desire to produce the following output data 



E0 



struct ures : 




of a hydrogen molecule. 




Fig. 3.2 - Two dimensional bounding box used for determining 
overlap percentage value. 
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* num_coords : number of coordinates generated for the two 

input planes. 

* new_coords (XYZ, num_coords) : coordinates generated by the surface 

construct ion process for the two~ planes. 

* new_conns(nurn_coords) : drawing instructions for each coordinate 

generated (SETPOINT, DRflWTO, DRAWPQINT). 

If the output data is in the form of triangular surface patches, an 
alternative data structure is required: 

* num_patches : number of surface patches generated for the 

input two planes. 

* new_CGGrds (XYZ) : new coordinates generated by the connection 

process. 

* patches (3, num_patches) : a 3 by num_patches array of triangles. 

B. THE ALGORITHM 

Our surface construction algorithm is composed of the following 
out 1 ined steps: 

( 1 ) Input and Inventory Compilation : 

The data structures defining the contours are processed to 
extract the pertinent data. This data includes the number of 

contours per plane, the coordinates defining those contours and 
the types of the contours. Additionally, two-dimensional 
bounding boxes are described about each contour for processing 
consideration in step 2. This compilation of data creates the 
data structures required for surface construct ion. 

(2) Overlap Deterrninat ion and Contour Item tapping : 

In this step of the algorithm, we determine which contours on 
adjacent planes have significant overlap, and which contours’ 
exteriors are near. This information is used to designate which 



contours should be connected via tr iangulat ions. The assignment 
of overlap is accomplished through the use of a value for the 
overlap percentage. This value is computed from the areas of 
the two-dimensional bounding boxes, as seen in Figure 3.2, of 
each contour. The overlap percentage is used to give priority 
to contour mappings that have the highest percentage of total 
overlap area. 

In this step of the algorithm, we also perform consistency 
checks for each contour pair. One such consistency check is 
executed using the contour interior specification and the 
overlap percentage value. Contour interior specifications are 
assigned as the value of a contour with respect to its immediate 
interior. As such, a contour is LOW valued if it is taken from 
the exterior of a solid object, such as the skin of an apple. 
Conversely, a contour is HIGH valued if its immediate interior 
is non-solid. Using these pieces of information, we are able to 
eliminate contour mappings of high overlap percentage which 
result in erroneous approximat ions of the original 
three-dimensional sol id. 

To illustrate the application of this consistency check, let us 
consider the mapping example for Figure 3.3. Here we are 
presented with a set of contours taken from a solid cone 
standing within a hollow cone. In this case, contour 1 on plane 
1 has a high overlap percentage with contour 2 on plane 2. 
However, since contour 2 on plane 2 is low valued with respect 



plane 1 
plane 2 




Fig. 3.3 - Example of consistency check using item interior 
specifications with overlap percentage values. 
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to its solid interior and contour 1 on plane 1 is high valued, 
this mapping can be eliminated. 

The interior specifications are also used to determine whether 
the mapping is interior to interior or exterior to exterior. fin 
int er ior to inter ior mapping is one which maps the 
interior of one contour to the interior of another contour. 
This form of mapping is indicative of contours taken from a 
surface with a shallow gradient, i. e. , a surface where the 
mapped contours are of similar size and shape, and where the 
contours have significant overlap. fin ex ter ior to 
ax ter ior mapping is one which maps the exterior of one 
contour to the exterior of another contour. This form of 
mapping is indicative of contours taken from a surface with a 
steep gradient, i. e. , a surface where mapped contours are of 
dissimilar size and shape, and where the contours overlap 
percentage is slight. Interior to interior mappings are more 
common. The exterior to exterior mapping is indicated for cases 
of two contours with a low percentage of overlap and differing 
interior speci f icat ions (HIGH: LOW, or vice versa). 
r orm the Coordinate tapping for each flapped Contour Pair : 

For each coordinate pair from step two, we form a complete 
coordinate to coordinate mapping. fi coord inate mapping 
is a tentative set of triangulation connections between the 
contour pairs. There are two procedures for determining this 
initial coordinate mapping. The procedure used is dependent on 
the type of mapping found for the paired contours in the 



previous step (interior to interior, or exterior to exterior). 
Additionally, both procedures try to form tr iangulat ion segments 
of shortest length, similar to the Christiansen algorithm. A 
general statement of this selection process is that we are 
trying to map coordinate i of contour n, plane 1 to coordinate j 
of contour rn, plane 2 such that the distance between the two 

coordinates is minimized. An additional qualification to this 
distance minimizing criterion is that coordinate connections do 

not cross, i. e. , coordinates 3 and 4 of plane 1 are not mapped 

to coordinates 6 and 5 of plane 2 respectively. 

(4) Continuity Recognition : 

The coordinate to coordinate mapping formed in step three is 
examined for continuity. Continuity, in this case, is defined as 
follows. ^irst, we form sets of coordinates from the coordinate 
mapping such that each coordinate of each set is constrained 
within a coordinate tolerance and within a distance range. The 

coordinate tolerance factor is a ratio of the number 

of coordinates in the larger contour divided by the number of 
coordinates in the smaller contour times a window value. The 
tolerance factor is used to group coordinates into a single set 
based upon their mapped coordinate number being within plus or 
minus tolerance of the last mapped coordinate added to the set. 
The tolerance sets formed are then compared for overlapping 
distance ranges. Any sets that have overlapping distance ranges 
are then merged. The merged set with the smallest distance in 



it is the set of coordinates for which connections should be 



generated. fill other coordinates are left unconnected. 

(5) Mapping Cancel 1 at ion ; 

Gnce we have decided to generate the connections for a part of a 
contour, we cancel any further mappings to that piece of the 
contour. This operation is required for partial mappings in 
which two or more contours on one plane are to be mapped to a 
single contour on another plane. This cancellation precludes 
connecting contour points which have already been selected for 
connect ion. 

(6) Connection Formation : 

We generate the coordinates for the triangulat ion connections 
specified in step four. "In between" coordinates, coordinates 
not directly mapped but within the tolerance factor for the 
connection mapping, are also added to the picture. The goal of 
the process is to form minimum area triangular surface patches 
for each segment of the mapped connection region. 

1 . Input and Inventory Compilation 

The input data to the algorithm consists of the contour 
descriptions for two adjacent planes of a three-dimensional solid. The 
purpose of this step of the algorithm is to segment this data into 
separate contour descriptions and to determine the individual 
characteristics of each contour. Figure 3.4 consists of two adjacent 
planes, each having three concentric rings of similar shaDe and 
continuity. Figure 3.6 consists of two closed loops on each of its 
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Fig. 3.4 - Example of multiple contours per plane on adjacent 
planes . 




Fig. 3.5 - Connection of Figure 3.4. 
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Fig. 3.6 - Example of a set of contours requiring partial mappings 
and an exterior to exterior mapping; (1,1) and (2,1) to (2,2). 




Fig. 3.7 - Connection of Figure 3.6, with contour interior values 
for each contour. 
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planes. Plane 1 has two small interior lobes, while plane 2 has one 
large surrounding contour with a small interior contour. The contour 
descriptions for these figures are composed of: 
the starting coordinate location, 
the total number of coordinates, 
the contour types, 
the interior values, and 

the contours’ two-dimensional bounding boxes. 

With the exception of the interior values, all of these character ist i cs 
are easily obtainable frorn the input data. 

The procedure necessary to obtain the contour interior 

speci f icat ions requires an evaluation of the data values lying along 
and interior to the contour (see Figure 3.3). If these values are not 
contained in the input data, a mechanism needs to be provided to allow 
for user specif icat ion of contour interior values. The range of 
interior values is HIGH, LOW or INDETERMINATE. Without this value the 
contour pairing operation encountered in the multiple contours per 
plane situation is difficult. In that case, some form of human 
interaction is necessary to designate which pairs of contours should be 
mapped together. If an interior value is not available, and the 
mapping situation is not complex, it can be set to INDETERMINATE 
without surface construction degradation. 

2. Overlap Deterrninat ion and Contour Mapping 

The overlap deterrninat ion and contour mapping procedure of tne 
surface construction algorithm is the process by which tentative 
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contour to 



contour mapping assignments are made. 



The contour 



characterist ics which are necessary for this procedure are the 
two-dimensional bounding boxes and the contour interior specifications. 
This mapping process is the key component in the disambiguation of 
multiply paired contours. 

The overlap determinat ion and contour mapping procedure is 
accomplished in the following manner. First, the two-dimensional 
bounding box of each contour on plane 1 is compared for overlap with 
the two-dimensional bounding box of each contour on plane £. The 
coordinates which define these bounding boxes are the minimum and 
maximum X and Y coordinates from each of the contour descript ions. 
(Additionally, these coordinates are adjusted by a constant value to 
promote overlap for exterior to exterior mapping situations.) From 
this operation, a table called the overlap table is produced. It is a 
two-dimensional table that contains a value for each possible pairing 
of contours between the two planes. The value recorded in each table 
entry indicates the extent to which each contour overlaps. If there is 
no bounding box overlap for a pair of contours, a value of 0.0 is 
recorded in the table. If there is overlap, the value recorded in the 
table represents the percentage of overlap with the larger of the two 
contours. This value is computed by dividing the area of the bounding 
box overlap by the area of the bounding box of the larger contour. 

After the overlap percentage has been computed for a contour 
pairing, it is used in conjunction with the interior specifications to 
determine the mapping type for the contour pair. An interior to 
interior mapping is indicated when a high percentage of overlap 
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(greater than 10’/.) exists for a pair of contours. A consistency check 
for matching interior specifications is performed for every pair of 
contours that exhibits this high an overlap. The consistency check 
requires that each contour pair have either HIGH:HIGH, LOW:LOW, or 
INDETERMINATEranything (HIGH or LOW) interiors. Contour pairings with 
high overlap but inconsistent interior speci f icat ions result in an 
adjustment to the overlap table of 0.0 percentage of overlap. An 
exterior to exterior mapping is indicated when the overlap percentage 
is low (less than 10’/.) and item interiors are non-matching. Finally, 
all contours with low overlap percentages and matching interiors are 
zeroed in the overlap table. 

Figures 3.8 and 3.9 graphically represent the overlap 
determ inat ion and contour mapping for Figures 3.4 and 3.6. Included in 
these figures are the overlap tables produced by this procedure. The 
table in Figure 3.8 shows three valid overlap percentages for three 
different contour pairs: (1,1) - (1,2), (2,1) - (2,2), and (3,1) - 

(3.2) . Four of the entries have been zeroed by the consistency check 
mechanism. Without this capability, high valued overlap percentages 
would appear in the overlap table with human interaction required for 
their disarnbiguat ion. The table in Figure 3.9 shows two high overlap 
percentages and two low overlap percentages. This data indicates that 
contours (1, 1) and (2,1) both map interior to interior with contour 

(1.2) . The low overlap percentages indicate that contours (1,1) and 
(2,1) map exterior to exterior with contour (2,2). 
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Fig. 3.8 - Bounding boxes and overlap table produced for Figure 3.4 
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Fig. 3.9 - Bounding boxes and overlap table produced for Figure 3.6 
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3. Form the Coordinate Mapping: Interior to Interior 



The coordinate mapping formation procedure for each coordinate 
pair having a non-zero overlap (in the overlap table) begins with the 
pair having the largest overlap percentage. All remaining steps in the 
surface construction algorithm are carried out on this pair before the 
next pair of contours is considered for mapping. Mapping paired 
contours is on a largest to smallest overlap percentage ordering. 
Since exterior to exterior mappings are indicated only in situations 
where the overlap percentage is low, they are considered for mapping 
only after all interior to interior mappings have been performed. This 
study follows that ordering and completes the description of the 
interior to interior mapping process before considering the separate 
process necessary for exterior to exterior mappings. 

The first operation performed on an interior to interior 
overlap pair is the deterrninat ion of which contour is interior to the 
other. This assignment is accomplished by comparing bounding box areas 
for the contour pair and designating the contour with the smaller area 
as interior. Once the interior contour assignment has been made, the 
center coordinate of that contour’s bounding box is computed. 

The knowledge of the center coordinate of the interior contour 
is used in the following manner. For each coordinate of the inner 
contour, we determine which coordinate of the outer contour is closest 
to a vector drawn from the center coordinate of the inner contour 
through the coordinate of the inner contour (see Figure 3.10). We add 
the qualification that the outer coordinate selected by this procedure 
must be farther from the center coordinate than the inner coordinate. 




Fig. 3.10 - Vector radiating from center coordinate through the 



interior coordinate towards the outer contour for tentative mapping 
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Fig. 3.11 - Example of a case where tentative mapping coordinates 
and associated distances vary greatly. 
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Also, the outer coordinate must be on the same side of the vector as 
the inner coordinate. The outer coordinates selected by this mapping 
process are recorded as the tentative coord inate map 

coord inat e for each inner coordinate. We also record the 
two-dimensional distance from each inner coordinate to its tentatively 
mapped outer coordinate. The resulting data structure contains the 

mapped outer coordinates with the distance to the inner coordinate to 

which it is mapped. 

The tentative connection map for Figure 3.4 is very good. Due 
to the similarity in size and shape of the mapped contour pairs, there 
is very little variation in the mapped distance values and the 
coordinates selected for mapping appear sequential. On the other hand, 
it can be seen in Figure 3.11, that large variations in distance values 
result from this tentative mapping process, and mapped outer 

coordinates appear with large gaps in the sequencing. This is due to 
the dissimilarity of the contour pair; the inner contour is relatively 
simple and much smaller than the convoluted outer contour. The 
procedure used to delineate a correct mapping from this tentative 
mapping is described below. 

a. Continuity Recognition 

The continuity recognition procedure uses the tentative 
connection map and associated distances for a pair of contours to 
determine the set of coordinate mappings that should be made for that 
pair. In the previous step of the algorithm, we produced the tentative 
connection map for all of the coordinates of the inner contour. This 
provides a rough approximation of the final mapping, but it must be 
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noted that not all of the inner coordinates need be involved in the 



final mapping for that pair. The continuity recognition procedure 
builds sets of coordinate mappings that are both continuous and of 
similar mapped distance range. These continuity sets are then used to 
determine the coordinate sequences that should comprise the final 
connection mapping. 

The first step in this procedure is to assign each 
coordinate pairing of the tentative connection map to an initial 
continuity set. This is accomplished by stepping through the 
coordinates of the inner contour in sequence and comparing each 
coordinates’ mapped outer coordinate to the last coordinate added to 
the last created continuity set. If that coordinate is within a 
tolerance factor of the last coordinate added, it is added to that set. 
If the coordinate in question is not within tolerance, a new set is 
created with that coordinate mapping as its start. The tolerance 
factor used is a ratio of the number of coordinates in the outer 
contour divided by the number of coordinates in the inner contour times 
a window value. (The window value is discussed in the next chapter. ) 

To illustrate this continuity set assignment, let us refer 
to the example in Figure 3.11. Here, the tolerance factor is 1® 
coordinates. The last coordinate considered is inner coordinate number 
24. The next coordinate considered is coordinate 25, which is mapped 
to outer coordinate 53. This coordinate is within the tolerance factor 
of 10 and is added to the last created continuity set. Inner 
coordinate number 26 is mapped to outer coordinate 69. This outer 
coordinate is outside of tolerance with the last coordinate added and 
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therefore, a new continuity set is created with this coordinate mapping 
as its start. 

This initial step of the continuity recognition process is 
a fast method for aggregating coordinate map pairs. In addition to 
building the initial continuity sets for the tentative mapping, we keep 
track of the minimum and maximum mapped distances for each continuity 
set. These values are used for merging continuity sets in the next 
step of the process. 

The initial sets generated for Figures 3.4 and 3.6 are of 
particular interest. This step of the continuity procedure placed all 
of the tentative mappings for the coordinate mapping pairs for Figure 
3.4 into a single set. This can be attributed once again to the 
contours’ similar shapes and sizes. On the other hand, coordinate 
mapping pairs for the mapping (1,1) - (1,2) of Figure 3.6 resulted in 5 
initial continuity sets with varying distance ranges (see Figure 3.12). 

□nee the initial continuity sets have been created for a 
contour pairing, we merge any sets that have overlapping mapped 
distance ranges. This merge process reduces the total number of sets 
and further aggregates the coordinate pair mappings to sets with 
coordinate number continuity and distance range similarity. In 
reference to our examples, no continuity set merge was required for 
Figure 3.4 due to its singular initial continuity set. Figure 3.12 
shows the initial sets with distance ranges and the merged sets with 
distance ranges for the contour pairing (1,1) - (1,2) of Figure 3.6. 
In that figure, the 5 initial continuity sets have been merged into 3 
sets of non-over 1 appi ng distance range. 
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Fig. 3.13 - Bounding box overlap for exterior to exterior mapping. 
Only the coordinates within the overlap area are mapped. 
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After we have merged continuity sets, we need to determine 
which of those sets of coordinate mappings is the one that should be 
used for connection formation. The choice is clearly" the set with the 
smallest distance range. With this decision, we validate all 
coordinate pairings that are members of this smallest distance set, and 
cancel all other coordinate pairings for that set of contours, 
b. Mapping Cancellation 

The validated coordinate connection map for the contour 
pair has significance beyond indicating which coordinates need to have 
connection segments generated. It also indicates "filled" connection 
positions. By filled we mean that once we have formed connections to a 
coordinate segment of a contour, that segment should not be reused for 
any further mapping that occurs for the two current, adjacent planes. 
This mapping is both checked and recorded at this stage of the 
algorithm. Mapping cancellation examines the coordinate mappings for 
which a validated mapping has been assigned. If either of the two 
coordinates, inner or outer, has been assigned to a higher priority 
mapping for this pair of planes, then that mapping is cancelled. Once 
these connections have been struck from the connection map, all 
remaining validated connections are recorded as filled. 

An additional tasking of this cancellation process concerns 
whether the mapping of either contour resulted in all coordinates 
defining that contour being included in the mapping. In that case, all 
other possible pairings with the completely mapped contour are 

cancelled. This is accomplished by zeroing the overlap on that 
contour’ s row or column of the overlap table. 
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c. Connection Formation 

When the above steps have been completed for a pair of 
contours, the remaining process of generating the* appropriate line 
segments is relatively simple. The final coordinate mapping for the 
inner contour is examined for continuous segments of validated 
connections. When a continuous segment is defined, the beginning and 
ending coordinates of that segment (for both the inner and outer 
contours) are used as boundary pointers for connection formation. The 
coordinates in between those pointers are stepped through one at a time 
by a process whose purpose is to generate the minimum area triangular 
surface patch, as defined in our introduct ion. The surface patch is 
formed by using a line segment from one contour as the triangle’s base, 
and a coordinate from the other contour for the triangle’s third point. 
The minimum area selection is accomplished by a procedure that chooses 
the next line segment between the contours that is both the shortest 
and within the mapping specified for the two contours. This is 
identical to the heuristic used by Christiansen in [CHRIS]. Differing 
coordinate rates between the two contours are taken care of by using 
the coordinate ratio (from the continuity tolerance factor) between the 
contours. This ratio allows the process to generate several line 
segments emanating from a single coordinate when there is a coordinate 
rate differential between two mapped contours. The lines generated by 
this procedure for Figures 3.4 and 3. £ are shown in Figures 3. Z ard 






3. 7, respect i vely. 



4. Form the Coordinate Mapping: Exterior to Exterior 

We begin the exterior to exterior mapping process at the same 
point of the algorithm where we departed in the description of the 
interior to interior mapping process. In keeping with our ordering 
criteria for mapping contour pairs, we examine the contour pa: - 
requiring an exterior to exterior mapping which has the highest overlap 
percentage in the overlap table. fill remaining steps of the algorithm 
are carried out on this pair before the next pair of exterior to 
exterior contours, in largest to smallest overlap area, is considered. 

In Figure 3. 13, we are presented with an enlarged view of the 
bounding box overlap area of the contour pairing (1,1) - (S,£) of 
Figure 3.6. This area of overlap contains all of the coordinates from 
both contours which are involved in the connection mapping. ~~he first 
operation performed on an exterior to exterior mapped overlap pair is 
the determination of the set of coordinates in both contours that ;= 
within the overlap area. The contour with the smaller number o* 
coordinates in the overlap area is used in the formation of a 
connection mapping between the contour with the larger number of 
coordinates in the overlap area. The basis for this connection map is 
the det erminat ion for each coordinate (in the smaller coordinate set 
contour) of the coordinate in the other contour coordinate set that is 
the shortest distance away. This det erminat ion is a simpler version of 
the distance minimizing process for connection set assignment c* 
interior to interior mappings. The product of this process is the 
connection map for the pair of contours. The use of continuity sets 
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is not necessary for exterior to exterior mappings due to the 
relatively srnall number of coordinates which comprise the connection 
set. 

Once we have generated this connection set, we use the same 
mapping cancellation and connection formation procedures as described 
for the interior to interior mappings. The connection formation 
procedure again uses the connection set mapping to find continuous 
segments of validated coordinate assi gnrnent s. The continuous segment 
thus defined is used to form triangular surface patches for all lire 
segments and coordinates within that segment. The final connection 
formation for the exterior to exterior mappings, (1,1) - (2,25 and 
(2,1) - (2,2) of Figure 3.6, are shown in Figure 3.7. 
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IV. ALGORITHM HEURISTICS PND LIMITATIONS 



In the preeeeding chapter, we presented an explanation of our 
algorithm for surface construct ion. Particular attention was devoted 
to our algorithm’s handling of the multiple contours per plane and 
partial contour mapping problems. It must be emphasized, however, that 
our algorithm does not provide a complete solution for all sets of 
contour surface data. In this chapter, we investigate some of the 
limitations of our algorithm. In order to do that, we must first 
discuss the heuristics employed by that algorithm. 

A. HEURISTICS 

Our algorithm utilizes three heuristics which are essential for the 
correct connection of planar contours. These heuristics were presented 
briefly in the last chapter, but we feel it is necessary to explain 
more fully their application and interaction regarding the contour 
mapping problem. 

1. Overlap Percentage Minimum 

In step two of our algorithm, we determine the percentage of 
overlap between contours on adjacent planes. These percentages are 
then considered in a consistency check for matching contour interior 
speci f icat ions. The heuristic in question, the -overlap percentage 
minimum, is applied in the final phase of this contour pairing 
procedure. Contour pairs having an overlap percentage value above the 
overlap percentage minimum, with matching interior specif icat ions, are 
designated for interior to interior mapping. Contour pairs having 
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non-zero percentages below the overlap percentage minimum, with 
non-rnatching interior sped f icat ions, are designated for exterior to 
exterior mapping. fill other contour pairs are disregarded. 

The value we have utilized for the overlap percentage minimum 
is ten percent. We found, through exper irnentat ion, that the assignment 
of this value resulted in the greatest number of correct contour 
pairings. Some contour pairs which should be mapped, however, are 
disregarded for mapping because of this selection (of 10’/-) for the 
overlap percentage minimum. In Figure 4.1, we are presented with an 
example of such a situation. In that figure, we have a pair of 
contours with matching interior specifications (HIGHsHIGH) , and having 
an overlap percentage less than ten percent. By our heuristic, this 
contour pair would not be considered for mapping, and would remain 
unconnected. 

One possible solution to this problem would be a mechanism 
which used a relaxation procedure to force a mapping between the pair 
of contours. This mechanism could be selected by the user to designate 
contour pairs for mapping which would otherwise be disregarded. If 
applied to the mapping situation of Figure 4.1, an appropriate 

connection could be generated. 

2. Boundary Tolerance Percentage 

The next heuristic to be discussed comes into play in the 
initial two steps of our algorithm. Specifically, the two operations 
involved are the determinat ion of contour item two-dimensional bounding 
box values, and the usage of those values for overlap determinat ion. 
As previously discussed, exterior to exterior contour mappings are 
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Percentage of overlap area < 10 % 



Fig. 4.1 - Example of a contour pair ’which should be mapped, but 
would be disregarded due to overlap percentage below the minimum. 




Fig. 4.2 - Example of contours 5 2D bounding boxes created strictly 
from the min and max X and Y coordinates. Resulting overlap = 0. 
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indicated for pairs of contours with a low percentage of overlap and 
non-matching interior sped f icat ions. In the initial development of 

our algorithm, we utilized the minimum and maximum X "and Y coordinates 
of the contour to describe its bounding box. We found, however, that 
in the majority of cases, these values resulted in zero percentage of 
overlap between contours which should be mapped. An example of this 
limiting of bounding box values can be seen in Figure 4.2. In that 
figure, we are presented with the contour pair from Figure 3. 13. In 
this example, it can be seen that limiting the bounding boxes for these 
two contours to their respective minimum and maximum X and Y coord irate 
values results in zero percentage of overlap. This is an 
unsat isfactory situation since the contours should be mapped. 

To remedy this situation, we adjust the bounding box values by 
a percentage to promote mappings in situations similar to that 
Figure 4,2. Once again, we are presented with the opportunity to 

utilize a relaxation procedure, prompted via user intervent ion, for 
mapping situations not included by this heuristic. A mechanism could 
be provided allowing the user to designate the bounding boxes fon 
individual contours, and thereby force a mapping between the desired 
set of contours. 

3. Tolerance ^ult ipl ier 

In an interior to interior mapping situation, a tolerance 
factor is used for the determination of the initial continuity set 
assignments. This tolerance factor is a ratio of the number of 
coordinates in the outer contour divided by the number of coordinates 
in the inner contour times a window value. The window value is a 
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constant which we found necessary for the selection of appropriate 
mapping connections. We chose to utilize a tolerance factor in this 
step of our algorithm, as well as in the connection formation 
procedure, because it provides an inexpensive means for restricting the 
search space in the selection of mapping connections. 

B. LIMITATIONS 

In the preceeding chapter, we demonstrated the capabilities of our 
algorithm, with emphasis on its handling of the problems of multiple 
contours per plane and partial contour mappings. We have found, 
however, that there exist contour mapping situations which cannot be 
handled by our algorithm. 

The first mapping situation concerns simple branching of one 
contour on one plane to two or more contours on an adjacent plane (see 
Figure 2.3). In this situation, we found that the application of our 
algorithm produces an incomplete contour mapping due to missing data. 
One possible solution to this mapping problem is the inclusion of a 
procedure for creating an introduced node similar to that described ir 
the Christiansen [Ref. 2j paper. This special case procedure could be 
selected automatically, or initiated via user interact ion. 

The next limitation of our algorithm manifests itself in situations 
where highly convoluted contours, with extreme narrowings, are mapped 
interior to interior. The problem here is due to the interior to 
interior algorithm’s dependence on the overlap region bounding box’s 
center coordinate for the tentative coordinate mapping. r or the 
portion of the contour near the center coordinate, the tentative 
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coordinate is fairly good. For the portion of a contour on the other 
side of a narrowing, where the center coordinate is no longer central, 
the tentative mapping is erroneous. The problem - comes when the 
tentative mapping is so bad that the continuity recognition procedure 
fails, and contour segments are incorrectly left unconnected. 

The solution to this problem is fairly simple and within the 
purview of our algorithm. If the convoluted contour is segmented at 
the extreme narrowings, it is possible to treat each open segment of 
the original contour as a separate contour. Using the original 
algorithm, we can generate centers for each new contour, and hence 
coordinate mappings, which result in a more correct approx i mat icn of 
the original three-dimensional object. The only capability lacking 
from our present algorithm is a mechanism for partitioning the original 
convoluted contour. This mechanism could be either user specified or 
automatic. The user specified option is favored due to the 
computational expense involved for automatic contour segmentation. 

The next limitation also concerns interior to interior contour 
mappings. In situations where sections of a contour tend to be near 
parallel with the vector drawn from the center coordinate of the inner 
contour, erroneous mappings result. Pn example of this situation can 
be seen in Figure 4.3. For those segments of the outer contour which 
are nearly perpendicular to the tentative connection vector, an 
appropriate connection map is generated. Ps the contour segrent 
becomes more parallel to this vector, the tentative connections 
generated begin to falter. 
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Fig. 4.3 - Example of situation resulting in an erroneous tentative 
coordinate mapping where contour segment becomes near parallel 
with the tentative connection vector. 




Fig. 4.4 - Example of a situation where two contours are mapped 
interior to interior which would result in an incomplete mapping. 
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The remedy to this problem is very similar to that for the previous 



situation involving highly convoluted contours with extreme narrowings. 
Segmentation of the original contour into several open segments, which 
could be mapped separately, would greatly improve the quality of the 
tentative coordinate mapping. Once again, user intervention is the 
preferred method of contour segment at ion. 

The final problem situation to be discussed concerns interior to 
interior mappings where the inner contour is not contained in the outer 
contour. This situation would result from contour data taken from a 
torus, such as a doughnut. An example is illustrated in Figure 4.4. 
The problem with this mapping situation results from the use of the 
tentative connection vector emanating from the center of the inner 
contour. Since the center coordinate of the inner contour is displaced 
from the center coordinate of the outer contour, tentative mappings are 
generated only for that section of the outer contour which is on the 
same side of the tentative connection vector (see Figure 4.4), “"he net 
result is a partial mapping of two contours which should be totally 
connected. 

A practical solution to this mapping problem, which could be 
readily adapted to our algorithm, is described in the Christiansen 
paper [Ref. 23. In mapping situations wnere contours to be mapped are 
not mutually centered, Chr ist iansen recommends a translation procecure 
onto a unit square, centered at (0,0). The principle of this process 
is to translate the two contours in such a manner that they become 
mutually centered within the unit square. Application of the interior 
to interior algorithm at this point would result in the desired 



mappings. Tentative mappings would be generated for the contours’ 
original coordinates, thus allowing the appropriate connections to be 
formed in the final step of the algorithm. 

C. SUMMARY 

It has been the purpose of this chapter to investigate the 
limitations of our algorithm, and provide practical solutions where 
possible. Additionally, to a lesser degree, the heuristics employed by 
our algorithm have been explained to increase the understanding of t u e 
reader. We feel that our algorithm provides a solution to the multiple 
contours per plane and partial mapping problems, but must concede that 
it is not a total solution to the problem of surface construction "rom 



planar contour data. 



V. CONCLUSION 



It has been the goal of this paper to describe a new algorithm for 
the surface construct ion of a three-dimensional object from a set of 
that object’s planar contours. The greatest part of this paper has 
been devoted to the capabilities of our algorithm, specifically, its 
handling of the multiple contours per plane and partial contour mapping 
problems. We have included a discussion of the limitations encountered 
thus far by our algorithm for specific problem mapping situations. 

In view of the limitations presented, we must comment that our 
algorithm does not, in its present form, provide a complete solution to 
the contour mapping problem. Further development is required to 
alleviate the problem areas discussed in Chapter IV. It is probable, 
however, that the correction of these algorithmic shortcomings will not 
ensure a complete solution to the contour mapping problem. We foresee 
that in some situations either user interaction or an alternative 
approach may be required. 
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APPENDIX 



PSEUDO-CODE ALGORITHM DESCRIPTION 



A pseudo-code description of our surface construct ion algorithm 
fol lows: 

FACE_IT 

{ 

Input the coordinates for two adjacent planes. Make a local copy 
of the coordinates. 

DEL I NEATE_ INVENTORY 
{ 

Take inventory of the contours in the coordinate sets. This 
inventory determines the total number of contours for each plane 
and records where each contour begins and ends. 

> 

TYPE_ INVENTORY 
■C 

We determine the contour type of each contour in each plane. 
There are three possible types: CLOSED_LCCP, OPEN_SEGMENT, and 
SINGLE_POINT. 

> 

BQUND_ INVENTORY 
i 

Determine the rectangular, two-dimensional boundary of each 
contour. Increase those boundaries by a constant to increase 
the possibility of detecting appropriate exterior to exterior 
mappings. 

> 

INTERIOR JDETERMINATIQN 
{ 

Determine whether the interior of each contour is HIGH or LOW 
valued with respect to the current contour level. This value 
can be assigned i nteract 1 vely in cases where the information 
to make this determ i nat ion is not available. These values 
are used in a consistency check for selection of contour pairs 
for mapping. 

> 

QVERLAP_DETERM I NATION 
i 

Compute the overlap table for the contours of both planes. The 
values in the table are the percentage of overlap for each 
possible contour pair on the adjacent planes. If there is no 
overlap, a value of 0.0 is recorded. 



Contour mapping types are also assigned at this step of the 
algorithm. Contour pairs with a HIGH percentage of overlap, 
matching interior specif icat ions (HIGH:HIGH, LGW:LGW, or 
INDETERMINATE: anything) are assigned interior to "interior type 
mapping. Those pairs with a non-zero overlap percentage, below 
10’/., with non-matching interiors are assigned exterior to 
exterior mappings. All other contour pairings are zeroed. 



CONNECT I ON JJETERM I NAT I ON 

This step of the algorithm orders the pairs to be mapped, and 
forms connections for the assigned types of contour mappings. 
This step is detailed below. 

> 



> /* end of FACEIT */ 



CGNNECTIGN_DETERMINATIGN 

while .true. 

{ 

Find the largest overlap percentage in the overlap table. If the 
largest value = 0.0 then QUIT. 

If the contour mapping indicated by this largest overlap value is 
exterior to exterior .... 

EXTER IQR__TQ_EX TER I OR _MAPP I NG 

Determine the set of coordinates in each contour that are in the 
overlap area. 

For the contour of the overlap pair that has the least number of 
coordinates, find the minimum distanced coordinate of the other 
contour. 

Assign all coordinates within the overlap region to the correc- 
tion set. 

> /* end of EXTERIGR_TG_EXTERIGRjMAPPING */ 

else 

/* perform an interior to interior mapping */ 

INTERIQR_TG_I NTE R I0R_MAPP I NG 

< 

Determine which contour of the pair is interior. This assignment 
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is based upon which contours’ bounding box is smallest. 

Compute the center coordinate of the inner contour’s bounding 
box. Check to make sure that this point is inside the contour. 
If it is not, the contour needs to be partitioned*. 

For each coordinate of the inner contour, determine the coord- 
inate of the outer contour which is closest to a vector drawn 
from the center coordinate through the coordinate of the inner 
contour. Store the coordinate as the connection map coordinate 
for the inner contour. Also, record the mapped distance from 
each inner coordinate to its mapped outer coordinate. 

RECOGNIZE_CONTINUITY 

i 

/* Determine continuity sets in the two contours using the 
the connection map and associated distances. *7 

INITIAL_CONTINUITY_SETS 

Assign the coordinates of the connection map to a con- 
tinuity set based upon whether each consecutive coordinate 
is within a coordinate tolerance factor. ""his tolerance 
factor is a ratio of the number of coordinates in the 
outer contour divided by the number of coordinates in the 
inner contour. 

> 

I N I T I AL_SET_D I ST ANCE. RANGES 

Determine the minimum and maximum distance ranges for each 
of the continuity sets. 

> 

CONTINUITY_SET_MERGE 

•C 

Merge any continuity sets that have overlapping distance 
ranges, maintaining the distance range for any merged set. 

> 

CONNECT I ON_SET_ASS I GNMENT 
i 

Assign coordinate connections for the coordinates of "he 
merged continuity set that contains the smallest distance. 
All other continuity sets are left unconnected. 

> 

> /* end of RECOGN I Z E_CONT I NU I TY */ 

> /* end of INTERIOR TO INTERIOR MAPPING */ 
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M APP I NG_C ANCELL AT I ON 

< 

Examine the coordinate mappings for which a connection has been 
assigned. If either of the two coordinates, inner contour or 
outer contour, has been used in a previous, higher priority 
mapping for this pair of planes, that coordinate mapping is 
cancelled. Once these filled connections have been struck from 
the connection map, all remaining validated connections are 
recorded as filled. 



CONNECTION_FORMATION 

Generate the connections for the validated coordinate map. Thi 
is accomplished by stepping through the connection map and 
forming coordinate connections where indicated. In between 
coordinates, those not directly mapped but within the tolerance 
factor for the connection mapping, are also added to the pictur 
The goal of the connection process is to form minimum area’ d 
triangular surface patches. 

> 

> /* end while .true, of CONNECT I ON_DETERM I NAT I ON */ 

> /* end of CONNECTION DETERMINATION */ 
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